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Wavelet-induced renormalization group for the Landau-Ginzburg model 
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The scale hierarchy of wavelets provides a natural frame for renormalization. Expanding the order parameter 
of the Landau-Ginzburg/$^ model in a basis of compact orthonormal wavelets explicitly exhibits the coupling 
between scales that leads to non-trivial behavior. The locality properties of Daubechies' wavelets enable us to 
derive the qualitative renormalization flow of the Landau-Ginzburg model from Gaussian fluctuations in wavelet 
space. 



1. INTRODUCTION 

Daubechies' wavelets||^,^ are an orthonormal 
basis that exphcitly separates scales. The basis 
functions are generated from a single set of func- 
tions ipt{x) by dyadic dilatations and translations. 
Each basis function is expressed as 
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where n G Z is the scale, x' E £" the posi- 
tion on a grid £" = (2"Z)-° with spacing 2", 
and t = l,...,nt — 2^ — 1 determines how 
the D-dimensional wavelet is composed of one- 
dimensional functions. The resulting wavelet 
fields naturally live on lattices £" of a dyadic 
multigrid. 

It was proved only by Daubechies that a com- 
pact basis with these features (and good analyti- 
cal properties) exists. It appears a natural basis 
to perform renormalization in. In this contribu- 
tion, we demonstrate how a wavelet basis can be 
used to exhibit explicitly the interscale coupling 
introduced by the ^I*"* term and how an approxi- 
mation to the renormalization flow in the Landau- 
Ginzburg/<f>^ model can be derived from this. 

2. WAVELET EXPANSION OF THE $4 
THEORY 

We consider the Landau-Ginzburg/$* theory 
|3| with a one-component real order parameter 
S{x), X e M^, governed by the Hamiltonian 
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Monte Carlo Simulation - Wavelet Fluctuation Strengths 




Figure 1. Wavelet fluctuation strengths at differ- 
ent scales as a function of the inverse temperature 
/3, measured in a Monte Carlo simulation 



with coupling constants vq and uq. The fleld is 
expanded in a wavelet basis: 

^(^)-EE E sl-\x')4"\x'){x)+S,{3) 

n t x'GC" 

Since wavelets have vanishing flrst moments, we 
use the real number to represent the overall 
magnetization of the system, while the S^"\x') 

s 



represent fluctuations at different scales. Fig. 

2\ 



shows their strength / sl."'\x') 



(2) 



as measured 

in a Monte Carlo simulation. As they show very 
clearly the location of the phase transition, we 
will focus on modelling these quantities. 

We thus make the variational ansatz that fluc- 
tuations are Gaussian and diagonal in wavelet 
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space: 

(4"^^(^'i)4"^^(4)> ='5„„„,<5,,,,,<5.i,., ^(^^^ (4) 

Note that local fluctuations in wavelet space still 
provide for a nontrivial correlator in position 
space, as the decay of the correlator is encoded 
in the relative strength of fluctuations at differ- 
ent scales n. 

The magnitude of the fluctuations ylj"'' is then 
found by the variational method from minimizing 
the free energy J- = U — S / (3 where the internal 
energy U is the expectation value of the Hamilto- 
nian H in the Gaussian ensemble (Q), and S the 
entropy of this ensemble. 

When calculating the internal energy, the self- 
similarity of the wavelet basis comes into play: As 
all basis functions are built from a single mother 
wavelet, the matrix elements of the Laplace op- 
erator have a simple scaling form: 

= 2-2«C*„*, . (5) 

Similarly, the four-point overlap integral occur- 
ring from the term has an approximate scaling 
representation. In this approximation, the effec- 
tive internal energy per site reads simply 

^ nt 

+ 'is^ + f 5^ (6) 
(5 — (Sq)) with the fluctuation sum 
^ = ^^2-^4") . (7) 

n t 

while the Gaussian terms are linear in the fluctua- 
tion strengths, the term introduces a coupling 
between scales. Since wavelets have compact sup- 
port, the overlap integrals and thus the nonlinear 
contribution is finite, as opposed to the situtation 
with a Fourier basis. The canonical dimension 2 
of the Laplace operator enters in the weight factor 
of the first term. 



Numerical minization - Wavelet Fluctuation strengths 
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Figure 2. Wavelet fluctuation strengths at differ- 
ent scales as a function of the inverse temperature 
/3, calculated from minimizing the effective inter- 
nal energy. 



Minimizing the free energy now with respect to 
the magnetization S = (So) yields 

5 = or 5= J-^-3A . (8) 
V 2mo 

Thus spontaneous symmetry breaking occurs 
when the fluctuation sum A exceeds the criti- 
cal value —ro/GuQ. Similarly, minimizing with 
respect to the ^^"^ results in 



f3 i 2-2"(-CtO + iro + 3uo{A + S^) ' ' 

This is an implicit equation for the fluctuation 
sum A that can be solved numerically. Fig. ^ 
shows such a solution. 

For uq — 0, (||) would be the wavelet repre- 
sentation of a exponentially decaying correlator 
with inverse correlation length y/rQ. We can see 
that the interaction leads to a redefinition of the 
inverse correlation length depending on the quan- 
tity A which is defined implicitly by this equation. 
In particular, the correlation length diverges at 
the same point as the second solution in (H) be- 
comes real, signifying the phase transition. 

3. RENORMALIZATION FLOW 

To derive the renormalization group flow of the 
theory, we apply the idea of removing the fine 
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Renormalization flow - 2D Landau-Ginzburg 
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Figure 3. Renormal- 
ization flow of the ef- 
fective Wavelet Haniil- 
tonian in two dimen- 
sions. Shown is the 
ro > 0; Wo < quad- 
rant of the coupling- 
constant plane. The 
arrows indicate the di- 
rection of the renor- 
malization map. The 
Gaussian fixed point is 
at the lower right cor- 
ner, the second fixed 
point at the center. To 
the right would be the 
symmetric phase, to 
the left the symmetry- 
broken phase. 



degrees of freedom to the minimization of the free 
energy, i.e., we perform the minimization scale by 
scale. 

Assume there exists a lowest scale n = 0, e.g. by 
already having eliminated all scales finer than n. 
We can then minimize the free energy with re- 
spect to keeping all coarser scales as vari- 
ables. Af^ will then become a function of the 
^j"'', n > 0. By reinserting this function in the 
expression for the free energy and expanding in 
a Taylor series, one reaches a new effective free 
energy, now only depending on n > 0, in 

which the new terms arising from the expansion 
can be absorbed (to some order) in a redefinition 
of the coupling constants. This yields a renor- 
malization flow as a mapping (ro,Mo) (rQ,Mg) 
in the coupling constant plane. It turns out that 
there exists the Gaussian fixed point ro = mq = 0, 
and in £> < 4 a second fixed point, corresponding 
to the Wilson-Fisher fixed point. 

The actual position of the fixed points still de- 
pends on the matrix element of the Laplacian 
^Ctt and thus on the type of wavelet. This 
shows a limitation of our approximation as we 
disregarded the correlation between neighboring 



wavelets which is basically determined by the ex- 
tent of the wavelet. 

4. CONCLUSIONS 

A wavelet expansion can be used to derive the 
properties of the Landau-Ginzburg model and its 
nontrivial renormalization flow even in a rather 
simple approximation. The crucial features we 
have made use of are scaling and self-similarity of 
the basis and locality of the basis functions. They 
enabled us to focus on the fluctuation strengths at 
different scales as the quantities of interest that 
govern the phase transition. The effective free 
energy of the system exhibits in a minimal way 
the coupling between different scales. 
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